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^ • Abstract 

^ ■ We study a one-dimensional large- U (N) gauge theory on a circle as a 

CN . toy model of higher dimensional Yang-Mills theories at finite temperature. 

To investigate the profile of the thermodynamical potential in this model, 
we evaluate a stochastic time evolution of several states, and find that an 
^ ■ unstable confinement phase at high temperature does not decay to a stable 

deconfinement phase directly. Before it reaches the deconfinement phase, 
it develops to several intermediate states. These states are characterised 
by the expectation values of the Polyakov loop operators, which wind the 
temporal circle different times. We reveal that these intermediate states are 
the saddle point solutions of the theory, and similar solutions exist in a wide 
class of SU{N) and U{N) gauge theories on including QCD and pure 
Yang-Mills theories in various dimensions. We also consider a Kaluza-Klein 
gravity, which is the gravity dual of the one-dimensional gauge theory on a 
spatial S^, and show that these solutions may be related to multi black holes 
localised on the S^. Then we present a connection between the stochastic 
time evolution of the gauge theory and the dynamical decay process of a 
black string though the Gregory-Laflamme instability. 



1 Introduction 



Understanding of the thermodynamics of gauge theories is a significant sub- 
ject in theoretical physics [1] . In the case of the pure Yang-Mills theory, 
due to the developments of the lattice gauge theory [2] and the holography 
in string theory [31 SI El El |7], its several properties have been confirmed. 
According to these studies, the theory is confined at low temperature and 
deconfined at high temperature, and a phase transition happens between 
these two phases. 

However we have not understood the natures sufficiently. For example, 
although it has been revealed that the deconfinement phase is stable at high 
temperature, it is unclear whether all unstable and meta-stable states, like 
domain walls and plasma balls [8], have been derived there. Besides, if such 
states exist, the profile of the potential, which connects these states, should 
be investigated to understand the dynamics such as the decay of these states 
to the stable deconfinement phase. 

In this article, in order to investigate these problems, we employ a "stochas- 
tic time evolution" , which is a kind of Monte Carlo calculation introduced in 
Sec. 14.11 We will apply this method to a finite temperature one- dimensional 
gauge theory at large A^, which we will explain in the next paragraph, as a 
toy model for higher dimensional Yang- Mills theories. Then we evaluate the 
stochastic time evolution of unstable states, e.g. a deconfinement configura- 
tion at low temperature and a confinement configuration at high temperature. 
Interestingly we find that the unstable confinement phase at high tempera- 
ture does not decay to a stable deconfinement phase directly. It develops to 
several intermediate states before it reaches the deconfinement phase. This 
evolution reflects certain characteristic aspects of the thermodynamical po- 
tential of this model. We will reveal that these intermediate states are related 
to the saddle point solutions of the potential, and similar solutions exist in a 
wide class of SU (N) and U (N) gauge theories on a circle including QCD 
and pure Yang-Mills theories in various dimensions. Thus these saddle point 
solutions may characterise the thermodynamical properties of these gauge 
theories too. 

Now we introduce the above one-dimensional U{N) gauge theory, which 
we mainly study in this article. The action is defined by 

S= /Vr('f]i(ArO'-ET[^''^']' 

\i=i i,j 
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Here is adjoint scalar (J = 1, ■ ■ ■ , -D) and Dt = dt — i[At, ■]. (5 and At 
denote the inverse temperature 1/T and the gauge field on the Euclidean 
time circle respectively. This model is a Kaluza-Klein reduction of a (D + 1)- 
dimensional pure Yang-Mills theory on S*^ x for small |9l |10] . 

We study this model at large [11]. We have three motivations to 
consider the large- limit. Firstly a large- iV phase transition occurs in this 
model, which is akin to the confinement /deconfinement phase transition in 
higher dimensional pure Yang-Mills theories [lOl [121 ESI [H] • Hence this model 
would be regarded as a toy model for the Yang-Mills theories. Secondly we 
can consider the dual gravity at large N [121 [13 [IS] . Lastly, since the large- TV 
limit is a kind of semi-classical limit, the role of the saddle point solutions 
would be emphasised. 

As we mentioned, we investigate this model by using the stochastic time 
evolution, and find the new saddle point solutions. In addition, we will also 
analyse them by using a, 1/ D expansion [TH [17] and reveal several relevant 
properties of these solutions including the critical temperatures, below which 
the solutions cannot exist. Besides we will compare the observables of these 
solutions from the Monte Carlo calculation and the 1/D expansion, and 
confirm good agreements. 

We will also study an application to the gauge/gravity correspondence 
[T2l [151 [16]. If we regard the temporal circle of the model ([1]) as a spatial S"^, 
this model at D = 9 can be obtained as an effective theory of A^ Dl branes 
winding this spatial at high temperature. Then we can construct the dual 
geometry by using these Dl branes |12j^ . Since the geometry involves the 
spatial S*^, the gravity becomes a Kaluza-Klein theory. We will argue that 
the saddle point solutions in the gauge theory would correspond to multi 
black branes localised on the S^. Indeed the existence of the dual phases in 
the gauge theories corresponding to these multi black branes in the Kaluza- 
Klein theories has been predicted in Ref. [16], and we will provide evidence 
for this conjecture. 

In addition, it has been shown in Refs. [TH [191 [2D] that the multi black 
holes on (connected by thin black strings) appear during a decay process 
of a black string due to the Gregory-Laflamme instability [2T] . We will argue 
that this phenomenon may be related to the stochastic time evolution of the 
unstable confinement phase of the gauge theory ([T]@. 

^Correctly, we will take a T-duality along the S*^ and consider DO brane geometries 
instead of the Dl branes as we argue in Sec. [S] 

^Indeed one motivation of this article is to understand the black string decay through 
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Such a connection to the gravity indicates that the stochastic time evo- 
lution must capture relevant properties of the large- gauge theories on S^. 
Since such gauge theories on are studied in various contexts: finite temper- 
ature gauge theories [1], higher dimensional extensions of the standard model 
in the phenomenology of particle physics [231 1211 [25] and string theories, we 
believe that our results will illuminate some non-perturbative aspects of these 
theories, especially in the dynamical problems. Concrete applications will be 
discussed in Sec. |6l 

This article is organised as follows. Although we find the saddle point 
solutions through the stochastic time evolution in the model ([T|), we first de- 
rive these solutions in four dimensional Yang-Mills theory and QCD through 
a weak coupling analysis in section m This argument will provide a basic 
idea of the saddle point solutions. In section [3l we study the model ([1]) 
by using the 1/D expansion, and reveal the natures of the saddle points in 
non-perturbative regime. In section HI we investigate the stochastic time evo- 
lution of the model ([T|) and show that the saddle points indeed appear as the 
intermediate states between the unstable confinement phase and the stable 
deconfinement phase at high temperature. In section \5[ we argue the role 
of these solutions in the gauge/gravity correspondence. Section [6] contains 
the conclusions and discussions. In appendix \^ we show some details of the 
analysis of the model In appendix [Bl we explain the details of our Monte 
Carlo calculation. 

2 Saddle Point Solutions in Gauge Theories 

In this section, we study saddle point solutions in SU{N) and U{N) gauge 
theories at finite temperature. The apphcation to the gauge theories on a 
spatial is straightforwarcif]. 

To investigate the gauge theories, the static diagonal gauge [1] 

(Af)., = (z,j = l,--- ,iV) (2) 

is useful. Here At is the temporal component of the gauge field, and does 
not depend on the Euclidean time t and has a periodicity 27r. If the theory 

the gauge theory, which is related to a ongoing project [52]. 

•^We have to take care of the periodicities of the fermions along the circle. It is fixed 
anti-periodic along the temporal circle, whereas we have to choose it either periodic or 
anti-periodic along the spatial circle. 
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is SU{N), they satisfy a constraint X^^^ = (mod 27r), although it is 
irrelevant at large A^, which we will study in later sections. 

The phases of the gauge theories are characterised by the expectation 
value of the Polyakov loop operator 

P ^ _TrP exp l^z Adtj = ]^ E • 

In adjoint theories, the phase is confined if {\P\) = and is deconfined if 

In the static diagonal gauge ([2]), the configuration of the eigenvalue {afc} 
determines the value of the Polyakov loop operator (|3]). Since can 
be regarded as the positions of particles on S^, this particle problem 
governs the phase structures of the gauge theories. Then we can interpret 
the mechanism of the confinement and deconfinement as follows. If repulsive 
forces between the eigenvalue {ofc} are dominant, the eigenvalues tend to 
spread on the and the stable configuration would be 

ttfc = 27rk/N + c (mod 27r), (4) 

and their permutations. Here c is a /c-independent constant, which is discre- 
tised as 2Tin/N {n G Zn) in the SU{N) case. Then \P\ = is satisfied and 
the confinement phase is realised. 

Oppositely, if attractive forces are dominant, the eigenvalues tend to 
clump and the configuration 

Q!fc = c (mod 27r) (5) 

would be stable. Sometimes quantum effects disturb this configuration and 
the eigenvalue distribution is smeared around c as we will see in Sec. [31 In this 
configuration, \P\ ^ is satisfied, and the deconfinement phase is realised. 
Therefore the phases in the gauge theories may be related to the forces for 
the eigenvalue {ak}- 

Then we intuitively notice that, if the attractive forces are dominant at 
certain temperature, the configurations with multiple mobs of the eigenvalues 
may be possible as a solution of the theory, if the attractive forces between 
the mobs are balanced. We call such a solution "multi-peak solution" and 
will find these solutions in several modelfl 

^Thc multi peaks of the eigenvalues of At are characterised by the expectation values 
of the generalised Polyakov loop (fT2|) too. 
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Figure 1: The effective potential for SU{3). We fix ^3 = —ai — a2 (mod 
27r) and plot ([6]) as a function of ai and 02 • The left is the three dimensional 
plot and the right is its contour line. In the right, ■, A and • denote the 
minimum, maximum and saddle point of the potential, respectively. 



2.1 Four Dimensional Pure Yang-Mills Theory 

We consider four dimensional pure Yang-Mills theory at high temperature. If 
we assume {ai} does not depend on the spatial positions, we obtain the one- 
loop effective potential for which would be valid at high temperature 

t»aqcd m, 

24 ^ VVtT 7r/mod2 / 

This potential implies that the attractive forces between the eigenvalues work, 
and indeed it is minimised by the configuration ([5]). Hence the system is 
deconfined at high temperature as is well known. 

We plot the potential for SU (3) in Fig. [1] Then we notice that this 
potential has several maxima and saddle points in addition to the vacuum 
solution The maxima are given by the confinement configuration 
The saddle points are given by {ak}k=i,2,3 = {tt + 27m/3, n + 27m/3, 27m/3} 
(n G Z3) and their permutations. At these saddle points, the eigenvalues are 
antipodal on the S^, and the attractive forces are balanced. 

Such multi-peak saddle points generally exist in SU (iV) and U (N) Yang- 
Mills theories if > 3. We show it by considering a two-peak solution as 



(6) 
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Figure 2: The two-peak solution of the potential (|6]). We put — 2M 
eigenvalues at and M eigenvalues at vr ± 5. The right plot is the potential 
for 5 (0 < 5 < vr). This plot displays that the two-peak solution (5 = 0) is 
meta-stable against this perturbation and is a saddle point of the potential 



an example. We arrange — 2M eigenvalues at and 2M eigenvalues at 
TT, and then this configuration becomes a solution of the action (|6]), since 
the attractive forces are balanced. Although this configuration is obviously 
unstable, it is meta-stable against certain Z2 symmetric perturbations. For 
instance, if we separate the 2M eigenvalues at vr to M eigenvalues at vr ± 5, 
then they tend to go back to vr for a small 5 as shown in Fig. |2l Thus this two- 
peak configuration is a saddle point of the effective potential (|6]). Similarly 
we can find various saddle point solutions by adjusting the positions of the 
eigenvalues. 

Particularly, if A^ is a multiple of an integer m, we can construct a 
symmetric saddle point solution, 

ai = 2Tik/m + c, {k - l)N/m< i < kN/m, {k = 1,2, ■ ■ ■ ,m). (7) 

We call such a Z„i symmetric solution "Z^ solution" in this articl^. Note 
that the uniform configuration is a special case of this solution with 
m = N but it is unstable against any perturbation and is a maximum rather 
than a saddle point of the potential 1^. 

^Similar Zm symmetric solutions may appear as stable states in gauge theories on 
coupled to adjoint fermions with the periodic boundary condition [26] . 
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TT 71 -71 7C -7C 7C 

(I) uniform (II) one-cut (III) non-uniform 

Figure 3: The distribution function p{a) in three states. 



2.2 Four Dimensional QCD 

We consider four dimensional QCD by adding A'^^ fundamental quarks to the 
pure Yang-Mills theory which we have argued in the previous subsection. 
In this case, the situation is slightly modified. The Nf fundamental quarks 
induce the one-loop potential for {ak} P, 
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N 



TT 



1) -1 

mod 2 



(8) 



in addition to the effective potential ([6]). This potential attracts the eigen- 
values at a/c = (mod 27r), and breaks the shift symmetry of ak- Even in 
this case, we can adjust the positions of the eigenvalues to obtain multi-peak 
saddle point solutions. 



2.3 Large- Gauge Theories 

We generalise the saddle point solutions in the previous subsections to large- 
gauge theories. At large A^, the difference between SU{N) and U{N) is 

irrelevant in our study and we ignore it. 

In the finite- A^ cases, the saddle points are characterised by the positions 

of the eigenvalues. At large A^, the eigenvalue distribution function [27] 

1 ^ 

i=l 
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-TT 71 -7C 7[ 

(IV) two-cut (V) gapless two-peak 



Figure 4: The distribution function p[a) for two-peak states with two gaps 
(=two cuts) (IV) and gapless (V). 



is a more convenient tool. Fig. |3] represents its profile, and (I) and (II) would 
correspond to the confinement (jlj) and the deconfinement ([5]) phases, respec- 
tively. (In (II), we turn on a quantum effect and the eigenvalues compose a 
cut@.) 

One additional ingredient of the large- gauge theories is the topology 
of the eigenvalue distribution function p{a). Although, both (II) and (III) 
in Fig. |3]have one peak, they have different topologies, i.e. (II) has one gap 
and (III) is gapless. Thus they must belong to different phases. Indeed a 
third order phase transition related to these two phases occurs in the two- 
dimensional lattice gauge theory [28| 129] . This type of transition is called 
"Gross- Witten-Wadia (GWW) transition" and occurs only at large A^. (The 
singularity at this transition point is resolved at finite A^.) In this way, we 
need to consider peaks, cuts and gaps in the eigenvalue distribution function 
in order to classify the phases. 

Now we consider multi-peak solutions in large-A^ gauge theories. The 
basic idea is similar to the finite- A^ case. If the deconfinement phase is stable 
owing to the sufficient attractive forces between the eigenvalues, the saddle 
point solutions may be possible by adjusting the positions. For example the 
Z2 symmetric two-cut solution like (IV) in Fig. H] may exist. In this solution, 
if the attractive forces are not so strong, the cuts would spread and the gaps 
might disappear. Then a GWW type transition might happen and a gap- 
less two-peak solution like (V) in Fig. H] might appear. In this way, various 

^We use the word "cut" for a region which is occupied by the eigenvalues and use the 
word "gap" for an unoccupied region in the distribution function p{a) according to the 
standard terminofogy of matrix model |27] . 
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topologies of multi-peak solutions would exist in the large- gauge theories. 

In this section, we have seen the multi-peak saddle point solutions in the 
gauge theories on S^. Similar calculations would be possible in weakly cou- 
pled SU (iV) and U (N) gauge theories, and these solutions may be generally 
founcjil. Note that related multi-peak solutions have been studied in several 
models by selecting appropriate setups so that these solutions become stable 
[26l [30l [3T] . What we have shown is that these solutions appear not only in 
these specific models but also in a wide class of gauge theories. 

Then a question is what the importance of these solutions is. We will 
study their roles in the following sections. 

3 Saddle Point Solutions in the One-dimensional 
Large- Gauge Theory 

In this section, we study the one dimensional gauge theory ([T]) at large A^. 
Although we have studied the multi-peak solutions in the gauge theories 
in the previous section, the analyses were just perturbative calculations in 
weak coupling. We investigate the non-perturbative aspects of the multi-peak 
solutions in the model ([T]) through a 1/D expansion. 

3.1 Phase Structure of the One-dimensional Gauge The- 



Before studying the multi-peak solutions in the one-dimensional large-A?^ 
gauge theory (Q, we summarise the phase structure of this model, which 
has been studied by using a. 1/D expansion 

We take a limit D, N ^ oo and (? — )■ with fixed A = g^ND, and then 
find a non-trivial vacuum, in which possesses a mass gap A, where A 
satisfies a relation 



^ Although the one-loop effective potential is described by the sum of the potentials 
between the two eigenvalues like ([6|), it generally involves higher order couplings like 
aiajak ■ ■ ■ in a strong coupling. Then the intuitive derivation of the saddle points based 
on the simple attractive force may not work. 



ory 




(10) 
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Thus we can integrate out and obtain an effective action [T^ : 



^°"<"-^)--g^^^ + f (11) 



DN^ 8A 2 ^ n\D 



n=l 



where is the generahsation of the Polyakov loop ([3]) which is 

defined by 



-.27r 1 



Mr 



Hence it satisfies 




p{a) = -\ l + \^u^e-^-^\. (13) 



The mass gap A is fixed such that this action is extremised. 

We evaluate this effective action in order to figure out the phase structure. 
Note that n„ for each n can be regarded as independent variables at large 
A^, as far as the eigenvalue distribution function p{q) does not have any gap. 
{{un\ have to satisfy a constraint p{a) > but it is irrelevant only if p(a) is 
gapless.) Then we obtain a solution for 



(n > 1), A = Ao = A^/^ (14) 



From eq. f|T3|) . we obtain p{a) = Xjlix and this solution corresponds to the 
uniform configuration Hence it is a confinement phase. 

Next we consider a low temperature regime, in which exp(— /3A) <^ 1 is 
assumed. There we can expand the action f lTT]) with respect to exp(— /3A) 
and obtain a saddle point for A as 

A = A^/3('l + ^f;e-"^^^^^Kp)+---. (15) 



n=l 



Then by inserting this solution to the action ffTTj) . we obtain an effective 
action at low temperature 

'^cff(^n.) _ 3 ~i/3 x-^ an I |2 , , ,4 , 
n=l 
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i-e-'^^^'^y h = ]-p\'/'e-'^'''\ (16) 



This expression implies that the uniform solution ( 1T4|) is stable if a\ > 0, 
which means that the temperature is lower than 

Tci = Y—^ + ---. (17) 
log -D 

Above Tci, fli becomes negative and the solution (HM becomes unstable. On 
the other hand, from Td, another branch of the solution 

ni = V-ai/26i, n„ = {n > 2), (18) 
Ai, ^ AV3 a _ ^^e^^'''') + ■■■ (19) 



appears, and this solution is stable above Td. From eq. f[T3|l . the eigenvalue 
distribution function becomes 



and this configuration is gapless and one-peak like (III) of Fig. [31 Note that 
this solution has broken the shift symmetry of a spontaneously and we have 
chosen the origin of a such that the peak of p{a) is at a = 0. In this way, a 
phase transition between this solution and the uniform solution (IM]) happens 
at Tci, and the order of this transition is second as shown in Ref. [H]. See 
Fig. E] too. 

The solution fl2Up cannot exist above 

Ai/3 

at which ui achieves 1/2 and p{a) becomes at a = ±7r. A GWW type 
third order phase transition occurs there and, above this temperature, the 
configuration possesses one gap like (II) of Fig. [31 There Un cannot be re- 
garded as independent variables anymore, and it is difficult to calculate m„ in 
this regime except T — Tc2 <^ A^^'^ as shown in Ref. [H]. On the other hand, 
if temperature is sufficiently high T ^ Tc2, At behaves as a free Gaussian 
matrix model with a mass M = \/2D/ /3Ai, where 

A = Ai = f 2Arj (22) 



11 



as we show the details in appendix |Xl The eigenvalue distribution becomes 




p(«) = — Re ./— . (23) 



Here the cut spreads on [— 2/M, 2/M], and it achieves the 5- functional dis- 
tribution ([5]) at T = oo. 

Therefore this model has at least three solutions: uniform f|T4|) . gapless 
one-peak fl2U]) and one-gap (= one-cut) fl2^ . and these solutions are stable 
in [0,Tci], [Tci,Tc2] and [Tc2,oo] respectively. We can also confirm these sta- 
bilities by evaluating the values of the effective action (fTTj) at these solutions, 

-^A^/^ (^T), (24) 



^^^AV3_|t, (T.<T<T.), (25) 

~^og(J-], (r,2«T). (26) 



(3DN^ 4 ^ VaV3 

Here the first and second actions are for the uniform and gapless one-peak 
solutions. (In this article, we put a prime to the quantities of non-uniform 
gapless solutions.) The third one, which is derived from ( I43p . is for the one- 
cut solution f l23|) at high temperature T ^ Tc2- Note that we can evaluate 
the effective action for the one-cut solution in T > Tc2 by using a Monte 
Carlo calculation of the effective action f lTTj) . and show that 5*1 is connected 
to 5*1/ at the GWW point Tc2- We sketch these relations in Fig. |5l From 
these three actions, we can read off the stabilities. In T > T^, 5*0 is higher 
than 5*1' or 5*1 and it indicates that the uniform solution is unstable in this 
region. 

Although the 1/D expansion predicts the two transitions at Td and Tc2, it 
has not been confirmed in the Monte Carlo calculation of the model ([1]). The 
Monte Carlo just indicates that some transition(s) occurs around a critical 
temperature ~ 1.30 at D = 2 and ~ 0.89 at D = 9 B (Here we have 



^ One reason for the difSculty of the decision of the critical temperature(s) in the Monte 
Carlo calculation is that Tc2 — Tci ~ 0{1/D) in the 1/D expansion and may be small. 
In addition, the transition is not sharp in the Monte Carlo owing to finite- iV effect. Note 
that there is a possibility that this transition is a single first-order transition rather than 
the two transitions [TH [T31 1311 [Ml- We will fix this problem in a future work. 
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5'eff(ltn)//3 

/9c2/3 /3c2/2 GWW(/^c2) uniform solution 




Figure 5: A sketch of the effective action Scs{un) of the several solutions in 
the one-dimensional gauge theory ([1]) through the 1/D expansion. The green 
and blue lines are for the gapless and the gapped solutions respectively. Note 
that this sketch is not rigidly plotted based on the equations, and the GWW 
point f3c2 is much close to /3ci in the actual plot at large D. 



used the unit A = g^N = 1, which we will use throughout the numerical 
calculation in this article.) However the details of the properties around 
might not be crucial in our study. 

3.2 Zm Symmetric Saddle Point Solutions 

Now we derive the multi-peak saddle point solutions from the effective action 
f fTTj) . At high temperature T ^ Tc2 , we can numerically derive an m-cut 
solution for a given set of the numbers of the eigenvalue {A^i, A^2, ■ ■ ■ , A^m} for 
each cut, if Ni ~ 0{N) and m <^ A^, as shown in appendix |Al In this section, 
instead of considering this generic situation, we study the Zm symmetric 
situation, which is a generalisation of the configuration ([7]). In this case, the 
analysis is much simpler and we can derive Zm solutions analytically. 

Recall that m„ can be treated as independent variables in the effective 
action f llip . Thus we have a classical solution for 

Un = 0, (ny^km, fc = l,2, ■■■)• (27) 
We substitute this solution into the action f llip and obtain an effective action 
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for A and Ukm (k = 1,2, 
DPP ~ m 



(28) 



This action is the same expression as the original action (11 ip by replacing 

Un — )• Unm and /3 — m/3, except the overall factor 1/m. Therefore if {Mi*^'"*} 
(n = 1,2, ■ ■ ■) is a solution of the action (II ip at an inverse temperature /3, 
Unm = u^n^ (n = 1 , 2, ■ ■ ■ ) is a solution of the new action (1251) at the inverse 



temperature m/3. 

By using this trick, we obtain the Zm symmetric gapless m-peak solution 
from (EUD, 

Um =^/-am/2bm, Un = (u ^ m) , (29) 
A =A^,^ = AV3 (^1 _ ^e-P''"'^ , (30) 

where has been defined in eq.( lT6|) . From eq. (|T3l) . the eigenvalue density 
becomes 



"("'^^ii+^r^"^'""')- '''' 

Indeed this configuration has a symmetry: a — )■ a + 27r/m. The above 
arguments imply that this solution exists in mTd < T < mTc2- At T = mTc2, 
a GWW transition happens and m gaps arise. Then this solution develops 
to the Zm symmetric m-cut solution, which is given by ( l44l) with A^; = N/m 
at high temperature T ^ mTc2- In this m-cut solution, Ukm 7^ and m„ = 
for n 7^ /cm (fc = 1, 2, ■ ■ ■ ) are satisfied. 

Similarly we obtain the observables in Zm solution. From eq. (!25|) . the 
value of the effective potential for the gapless m-peak solution (129!) at mTd < 
T < mTc2 is given by 



SzL 3r. 



^1/3 _ (32) 



2 



I3DN^ 8 Amh 
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At high temperature T ^ TnTc2, A and the effective action of the Zm sym- 
metric m-cut solution becomes 

from eq. fl2^ and fl26l) . Hence solution for larger m has a higher potential 
and is more unstable. See Fig. O We will compare these results with the 
Monte Carlo calculation in Sec. 14.31 

Lastly we argue the stability of Z^ solution. In the effective action (!28|l . 
Zjji solution is the stable solution for T > mTd. It indicates that Z^ solution 
is stable against the perturbation with respect to Unm- On the other hand, 
a„ (1 < 77, < m — 1) are negative at T > mTd, and u„ = (1 < < m — 1) 
is unstable in the original action f|TT]) . Thus Zm. solution is unstable against 
the perturbation with respect toM„(l<?T,<m — 1) and is a saddle point 
as we expected. 

4 Stochastic Time Evolution in the One-dimensional 
Gauge Theory 

We have shown that the model ([I]) has many saddle point solutions. However 
we have still not fully revealed the natures of the model ([T]). One missing 
property is the profile of the effective potential for {un} (or {tti}), which we 
have seen in the SU (3) case in Fig. [TJ Such a profile is crucial to investigate 
the decay process of an unstable state to a stable state. Although our effective 
potential is thermal and is not directly related to the actual real-time decay 
process, if the decay happens sufficiently slowly or adiabatically, the decay 
process may reflect the profile of the potential. 

Of course we cannot draw the potential like the SU{3) case, since we have 
infinite variable at large A^. However we possibly read off the relevant part 
of the profile of the potential through the following procedure. Suppose that 
we can integrate out adjoint scalar in the model ([T]) and obtain an effective 
action for the gauge field S'cfr(«i). Then we consider an unstable solution in 
the action 5'cfr(aj). (It corresponds to ▲ in Fig. [1]) We smoothly deform 
{ai} from this solution such that S'efr(tti) is becoming smaller. By repeating 
this deformation, {oi} may settle down to a stable configuration finally. (It 
may correspond to ■ in Fig. [TJ) Then we can speculate the profile of the 
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potential between the unstable solution and the stable configuration from 
the history of the deformation of 

We investigate this process by using a stochastic time evolution of a Monte 
Carlo calculation, which we will explain in the next subsection. Then we will 
see that this method captures several characteristic properties of the model 
([T]), which are related to the dual gravity as we will argue in Sec. O 

4.1 Overview of the Stochastic Time Evolution 

We introduce the stochastic time evolution, which is designed such that the 
above process is realised. We assign a discrete "stochastic time" s for and 
, which is distinguished from the Euclidean time t in the model ([T]). We 
take an appropriate initial configuration at s = 0, and, to gain the time s, 
we update ai{s) and Y\s) through the following rule: 

1. Set a trial configuration aj,triai(s + 1) = o^iis) + r, where r is a small 
random number. 

2. If S[ai^tTia.i{s + 1)] < 5'[aj(s)], we accept this trial configuration as 
ai{s + 1) = ai,trial(s + 1). 

3. Even if S[ai^tTia.i{s + 1)] > ^[^^(s)], we still accept ai,triai(s + 1) with the 
probability exp(— S'[ai^triai('S + 1)] + and, if it is rejected, we 
retain ai{s + 1) = ai{s). 

4. Update the scalar field Y\s) many times until they arrive at an equi- 
librium for the given configuration + 1)}, and use this state as 
Y'{s + 1). 

The last step is taken to focus on the dynamics of {tti} in Scsio^i)- Through 
this evolution, {ai} is deformed gradually such that the action Sesica) tends 
to be smaller as we intended. Note that this procedure for {at} is based on 
the Metropolis algorithm. 

The technical details of this Monte Carlo calculation are explained in 
appendix [Bl 
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4.2 Decay Patterns in the Stochastic Time Evolution 



We investigate the stochastic time evolution of the unstable states of the 
model ([1]). At T > Tc[^, we take the unstable uniform solution f lT^ (confine- 
ment configuration) as the initial state, and evaluate the evolutions repeat- 
edly by changing the temperature and random number. Then we observe the 
following two evolution patterns depending on temperaturj^: 

Direct Decay In this pattern, the unstable uniform solution directly evolves 
to a more stable one-peak state (deconfinement configuration). See the left 
plot of Fig. [61 In this figure, we plot the stochastic time evolution of \un\ 
with respect to the stochastic time s. Through eq. (fT2!) . \un\ roughly indicates 
the number of the peaks of the eigenvalue density p{a). If m peaks exist in 
p(a), Um tends to have a stronger signal relative to other m„ {n ^ m). If all 
Uri are close to zero, it indicates an uniform distribution. Indeed, in Fig. [H] 
(LEFT), \un\ ~ around s = indicates the uniform state f lT^ and the 
signal > \u2\ > ■ ■ ■ in the late history indicates the one-peak state. This 
pattern is mainly observed at lower temperature Tc < T < c{D)Tc. Here 
c{D) is a constant, which seems to depend on D as c{D = 2) ~ 3.5 and 
c{D = 9) ~ 5.0, although the change of the decay pattern at c{D)Tc is not 
sharp. 

Cascade Decay In this pattern, the unstable uniform solution first evolves 
to a multi-peak state. Then the peaks attract each other, and two of them 
collide and merge into one peak. By repeating such collisions, the number 
of the peaks decreases one by one, and it finally achieves the one-peak state. 
This pattern is dominant at high temperature c{D)Tc < T. The m-peak 
states with a larger m tend to appear at higher temperature. 

One example is shown in the right plot of Fig. [6] and Fig. [71 In this 
example, a three-peak state as shown in Fig. [3 appears around s = 300. 
Correspondingly a strong signal appears at |n3| in Fig. [6l Then the two 
peaks merge into one and it becomes a two-peak state around s = 800, and 
finally they collapse to a one-peak state around s = 1300. In Fig. [Q the 
strong signal of \u2\ around s = 800 corresponds to the two-peak state but 

~ 1.3 at = 2 and ~ 0.89 at £> = 9 in the Monte Carlo calculation of the 
model ^ as we mentioned at the end of Sec. 13.11 

^°The movies for the stochastic evolutions are available on 
http : / / www2 . yukawa . kyoto-u . ac . j p/~azuina/ mult i_cut/ index . html 
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Figure 6: Two decay patterns of the stochastic time evolution of the uniform 
solution f|T4l) . s is the stochastic time, and we take D = 9 and = 60. 
(LEFT) The direct decay observed at T = 2.0(~ 2.2Tc). (RIGHT) The 
cascade decay at T = 6.0(~ 6.7Tc). 
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Figure 7: The stochastic evolution of the eigenvalue distribution in the right 
plot (the cascade decay) of Fig. El e*"'' {k = 1, - ■ ■ , TV) are plotted. 
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Figure 8: The decay of an unstable one-peak configuration to an uniform 
state at T = 0.84(~ 0.9T,) with = 60 and D = 9. In the right, we 
compare A with the \/D expansion for the uniform solution (fT4|) . 

the strong signal of \u^\ around s = 1000 is merely owing to the asymmetry 
of the eigenvalue distribution during the decay. 

In this way, the multi-peak states appear as the intermediate states in 
the stochastic decay process of the unstable uniform stat^. These results 
indicate that the multi-peak states lie between the uniform and the one-peak 
configuration in the potential valley of Scsic^i)- Note that, although the 
stochastic time evolution is powerful to capture these qualitative features of 
the potential, we have not established a proper method to investigate them 
quantitatively. 

Similarly we attempt the stochastic time evolution of an unstable one- 
peak configuration at low temperature T < Tc, but it merely evolves to the 
uniform configuration as shown in Fig. [Hi This would be consistent with the 
fact that the uniform solution f lT^ may be an unique solution at T < Tc 0. 

Note that we have to take care of the topology of the eigenvalue den- 

^^We observe the similar decay patterns in the Monte Carlo calculation of the effective 
action too but we omit to show the results here. In this calculation, we have to find 
A, at which the action (jf f p becomes the maximum rather than the minimum, since A is 
originally a complex variable |14| . 

One issue is around T^. As we mentioned in the footnote [51 the order of the phase 
transition of the model ([T]) has not been fixed. If the transition is the first-order, meta- 
stable solutions must exist around Tc [HJ [131 [SH [33] ■ 
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sity p(a) in the Monte Carlo calculation. If we smoothly deform p{a) at 
large A^, its topology never changes. The topology change must accompany 
some singularity such as the GWW transition. Thus if we take the uniform 
configuration as the initial state, any gap never appears during the evolution. 

On the other hand, the Monte Carlo calculation is done at finite N, and we 
cannot define the topology of p{a) rigidly. Therefore we have to speculate 
the topology at large from the result of the Monte Carlo calculation. 
For example, no gap should appear in the distribution function p{a) in the 
stochastic time evolution starting with the uniform configuration, if it were 
in large N. Hence we should interpret that the multi (one) peaks of p{a) 
in Fig. [7] are connected by tiny non-zero eigenvalue densities. Since such 
tiny densities cost little energy, they may be irrelevant in the dynamics. 
However they may be relevant in the gauge/gravity correspondence, since 
the singularity associated with the topology change may correspond to the 
naked singularity in the gravity [22l [33l [Ml [35] . 

4.3 Comparison of the Monte Carlo with the 1/D Ex- 
pansion about Zm Solution 

We have studied the decay patterns in the stochastic time evolutions through 
the Monte Carlo (MC) calculation. By applying this method, we evaluate the 
observables of Zm solution, and compare the results from the 1/D expansion. 

To obtain a Zm solution in the MC calculation, we take an appropriate Zm 
symmetric configuration as the initial state of the stochastic time evolution. 
We evaluate this evolution and find that a Zm solution may be realised as 
a long life stat^ at the initial stage of this process at High temperaturj^. 
The solution finally decays to the one-peak configuration, which is stable at 
T > Tc. See Fig. [9] for the Z2 case. From the arguments below eq.( l3T]) . Z2 
solution has a strong signal at \u2\, and |'U2n+i| should be zero. It indeed 
happens in Fig. M until around s = 1000. Thus we believe that this is the 
Z2 symmetric two-cut solution. Then we compare the observables of the MC 
and the 1/D expansion, and these are consistent as summarised in Table [TJ 

^■'The terminology "long life" would not be proper, since this state is not meta-stable. 

^'*Wc tried to observe any signal for the critical temperature for Zm solution around 
mTc- However, the life time of the Zm solution around mTc seemed too short and we have 
not found any clear evidence. 
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Figure 9: The stochastic time evolution of a Z2 symmetric state. We take 
D = 9, = 60 and T = 6.0(~ 6.7Tc). Z2 solution may be realised between 
s = 100 and s = 1000. At that period, the observables are also stable as 
shown in the right plot. We evaluate the averages of the observables there 
and compare them with the results from the 1/D expansion as summarised in 
Table [1] Note that A in the right plot remains almost constant until around 
s = 2200, although m„ start evolving around s = 1000. We will discuss this 
reason in appendix |X1 

















MC 


3.3625(5) 


2.850(2) 


2.602(2) 


36.71(2) 


19.14(4) 


13.46(5) 


1/D 


3.22 


2.71 


2.45 


40.5 


20.3 


13.5 



Table 1: The comparison between the 1/D expansion and Monte Carlo cal- 
culation of the model ([T]) for solution. A^, D and T are the same as 
the data in Fig. O We use the unit A = 1. We obtain Z3 solution by a 
similar fashion to Z2 solution. We evaluate A defined in ffTOl) and the in- 
ternal energy E = —j^XogZ [36]. The 1/D expansion predicts A^,^ and 
Ez^/DN^ = 3T/(4m) + 0{1/D) from eq.([33D. 
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5 Multi-cut Solutions in the Dual Gravity 



5.1 The Dual Supergravity of the One-dimensional Gauge 
Theory 

We have studied the multi-peak solutions in the gauge theory ([T]). Now 
we consider their dual geometries through the gauge/gravity correspondence 
[21 [15] . In this subsection, we review the dual gravity of the model ([1]) based 
on the arguments in Ref . [12] . Roughly speaking, the eigenvalue distribution 
of {oj} is related to the positions of black DO branes on the spatial circle in 
a Kaluza-Klein gravity. 

We consider IIB string theory on x 5"^^ x S}^, where we have defined 
the inverse temperature as (^2 and the period of the spatial circle as L. We fix 
the periodicity of the fermions along 5*^ as periodic. We put N Dl branes, 
which wind this Sj^^ x S*^, and then the effective theory of the branes is given 
by the two-dimensional maximally supersymmetric U (N) Yang-Mills theory 
on the x Si- 

Wm=^ f^dt /"Vxr + 5^ i(Z}^F^)'-i5^[F^F^]2 + fermions 

92 Jo Jo L I,J 

(34) 

In this theory, the temporal Kaluza-Klein (KK) modes and fermions are 
decoupled, if the temperature is sufficiently high <^ {L/\2Y^^, where 
A2 = 92^. Then the theory is reduced to the one- dimensional model ([T]) 
with D = 9 and A = X2/P2- Here we have to identify 

2d SYM on S},^ x Si Id model on S} ([ID 

^ + R 

(35) 







t, /3, 






A, 


4°' 




gY, 



where "(0)" in the two dimensional fields denote the zero modes with respect 
to the temporal KK expansion. 

Through the gauge/gravity correspondence, the gauge theory (134]) at large 
N can be described by the black Dl brane geometry in the IIB supergravity 

m 
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with a dilaton and Ramond-Ramond (RR) potential. This gravity descrip- 
tion is vahd at strong couphng X2P2 ^ 1 so that the a' corrections are 
suppressed. In addition to this condition, L should be larger than the ef- 
fective string length (32^^/ X^^^ at the horizon so that the excitations of the 
winding strings along S}^ are suppressed. However interesting phenomena 
will occur when L is smaller than this effective string length. Although the 
IIB supergravity does not work in this regime, we can avoid this problem by 
taking a T-duality on the S}^ and go to the IIA frame. Here the original IIB 
and IIA theory are related as 

IIB T-dual on Si IIA 

Si ^ Sh (L' = (2vr)VL) , 

Dl brane winding Si ^ DO brane localised on Si, 

Thus the period of the dual circle, which we have defined as L' = (27r)^/L, 
can be large for a small L, and the IIA supergravity description works if 
L ^ /Jg^^/Ag^^. Under this T-duality, the black Dl brane fl5B]) is mapped to 

^ = F-^ fdt' + F-^{^ + dx[' + u'dn^^ , (38) 

where x'l = x'l + L'. This geometry represents that the DO branes, which 
are the T-dual of the Dl branes as in (1371) . are uniformly distributed along 
Sl,, and is called "smeared black DO brane". Note that the topology of the 
horizon of this geometry is Sl, x S*^, and this geometry is a kind of a black 
string in Kaluza-Klein gravity. 

However this geometry is stable only for small L' fl2\ [T6| |37] . Around 
L' ~ X^^"^ / IS^^'^ , owing to the Gregory-Laflamme instability [21], the smeared 
black DO brane becomes unstable. A first order phase transition occurs and 
the stable geometry for larger L' is given by a black DO brane localised on 
a point on Sl,. In this geometry, the horizon is also localised on Sl, and its 
topology is S^. This geometry is called "localised black DO brane" and the 
metric is approximately derived in Ref. [T6] . 

Let us compare this gravity result with the one- dimensional gauge theory 

1 /2 

([T]). Note that the gravity is valid at low temperature (32^ I/X2 , whereas 
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the gauge theory description is vahd at high temperature ^ (-^7^2)^^^, 
and these two parameter regions are different. Nevertheless we can see that 
the stabihty of the gravity is quahtatively consistent with that of the gauge 
theory studied in Sec. 13.11 

The T-duality (l37|l maps the gauge field A^i on the Dl branes to the 
adjoint scalar Y^, which represents the positions of the DO branes on Sj^,. 
Hence the eigenvalue distribution of A^i may correspond to the DO brane 
distribution on Sj^,. Through the relation (!35l) . the uniform configuration (HM 
and the one-cut configuration like (l23|l may be identified with the smeared 
black DO brane (138|) and localised black DO brane respectively. Recall that 
the uniform configuration is stable for large L, which corresponds to a large 
(3 in the model ([T]), and the one-cut one is stable for small L. Through the 
relation L' = (27r)^/L, this stability is consistent with the IIA supergravityF^. 

5.2 Gravity Duals of the Multi-cut Solutions 

Now we consider the dual geometries of the m-cut solutions in the gauge 
theory ([T]). Through the correspondence between the eigenvalue distribution 
of A^i and the distribution of the DO branes on S}^,, the m-cut solutions would 
correspond to m black DO branes localised on S}^,. Such multi black branes 
have been known as the solutions of the supergravities, and are unstable due 
to the gravitational attractive forces between the branes. It is consistent 
with the instability of multi-cut solutions in the gauge theory. Indeed the 
existence of the dual phases in the gauge theories corresponding to these 
multi black brane solutions in the Kaluza-Klein gravities has been predicted 
in Ref. [16]. The saddle point solutions in our article provide evidence for 
this conjecture. 

5.3 Black String Decay 

Interestingly such multi black holes (connected by thin black strings) appear 
in a real-time decay process of a black string [TSl HSl [20]. As we mentioned, 

^^The order of the phase transition is different. In the gravity, the first order transition 
is predicted |12j . On the other hand, the gauge theory predicts the two higher order 
transitions at Td and Tc2 through the 1/D expansion |14j . This is not inconsistent, since 
the order of the transition generally depends on the couplings. Besides the Monte Carlo 
calculation of the model ([T]) has not fixed the order of the transition at Tc- 
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r-i 

Figure 10: A sketch of the decay process of the meta-stable black string in 
Refs. [iHl EH]. The last step is just a speculation. The inside of the apparent 
horizon is depicted as the black region. In Ref. [19], an interesting fractal 
structure, which is akin to a low viscosity fluid, is observed but we omit it 
in this sketch. 



the black string is unstable at large L', and a localised black hole is ther- 
modynamically favoured. The authors in Refs. [181 llH] examine the decay 
process of a meta-stable black strin^^ by solving the Einstein equation nu- 
merically. They found that several points of the horizon are growing and 
other parts are shrinking owing to the Gregory-Laflamme instability. As a 
result, a sequence of black holes joined by thin black string segments appear 
along Sj^,. See Fig. [101 However the numerical analysis did not work for 
sufficiently long time and the final state of this process has not been found. 
Since this configuration is unstable, it might evolve to a single localised black 
hole joined by an extremely thin black string. 

Remarkably this time evolution of the black string is similar to the stochas- 
tic time evolution of the unstable uniform configuration in the gauge theory. 
The multi-peak state appears as the intermediate state in both cases. This 
agreement may imply that the dynamical stability in the gravity may be 
related to the thermodynamical stability in the gauge theorjo. (Recall that 
the stochastic time evolution may refiect the profile of the thermodynamical 
effective potential of the gauge theory.) 

Note that, it is expected that the thin black string segment appearing 
in the above decay process would be pinched off through a quantum effect, 
when the size of the horizon achieves the Plank scale. Although it is difficult 
to investigate such quantum effect in the gravity, it may be possible in the 

^^In the model in Refs. [T51[Tn], the Gregory-Laflamme transition is first order and a 
meta-stable black string solution exist around the critical point. 

^^In gravity, the connection between the dynamical stability and thermodynamical sta- 
bility is argued in Ref. |3H]. 
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gauge theory, in which the quantum effect may correspond to the finite- iV 
effect. Hence it may be important to study the model ([T]) further to uncover 
this problem. 

However there are several discrepancies between the black string decay 
[T8l [T9] and our stochastic evolution in the model ([!]). We will discuss them 
in Sec. El 

6 Conclusions 

We have found the multi-peak saddle point solutions in various gauge theories 
on S^. The important condition for the existence of these solutions is the 
attractive forces between the eigenvalues, which typically realise in case the 
deconfinement phase is stable. Since many maximally supersymmetric Yang- 
Mills theories are deconfined at any finite temperature, these solutions may 
exist there. Then their gravity duals may also exist in the supergravities and 
they might belong to a new class of geometries. To show it, we are studying 
the supersymmetric gauge theories and will report on a future work. 

Through the stochastic time evolution of the one- dimensional gauge the- 
ory ([1]), we found that these multi-peak states may exist as the intermediate 
states between the unstable confinement phase and the stable deconfinement 
phase at high temperature. This result indicates that, even in the real-time 
decay process of the unstable confinement phase, some related intermediate 
states might appear. It is valuable to investigate it further but one difficulty 
is the definition of the multi-peak states in the real-time formalism, since the 
theory does not have the temporal circle anymore. 

A natural generalisation of this study is the stochastic evolution of an 
unstable confinement phase in higher dimensional gauge theories. In this 
case, at can depend on the spatial positions, and this may lead to much 
richer structures involving the domain walls and plasma balls [8]. 

We can also take the 5*^ as the spatial circle as we considered in Sec. |5l 
Then we do not have any difficulty about the definition of the multi-peak 
solutions in the real-time formalism. One interesting application is to the 
phenomenology of particle physics. In some models, the gauge field of the 
spatial 5"^, which is an extra-dimension of our world, is identified as the Higgs 
field [231 IMl |25] • If we consider the time evolution of such models in early 
universe, certain multi-peak configurations might appear and probably play 
some roles. 
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We have also studied the apphcation to the gauge/gravity correspondence 
by identifying the eigenvalue distribution of the spatial gauge field and the 
positions of the D branes on the spatial circle in the Kaluza-Klein gravity 
through the T-duality. There the multi-peak solutions in the gauge theory 
may correspond to the multi black branes localised on the circle. 

In this case, we found an interesting similarity between the decay pro- 
cess of the black string [HI [191 120] and the stochastic time evolution in the 
gauge theory. In relation to it, we have three challenges. First one is the 
investigation of the thermodynamical properties of the two dimensional su- 
per Yang-Mills theory ( I34p . which was partially done in Ref. [3H]- Especially 
the decision of the order of the phase transition is important. At the high 
temperature, the analysis of the one dimensional gauge theory ([1]) may be 
valid, and there the 1/ D expansion predicts that the phase transition at Td 
is the second order transition followed by the GWW transition at Tc2 [E]. 
However, at the low temperature, the gravity calculation is valid and it pre- 
dicts that the Gregory-Laflamme transition is first order [12]. Thus the order 
of the transition would depend on the parameters and it will be important 
to confirm it directly by investigating the model ( 134]) . In addition, the explo- 
ration of the relation between the cascade in the stochastic time evolution 
and the order of the phase transition may be also interesting. 

The second challenge is the study of the real-time evolution of the uniform 
configuration of A^i in the two dimensional super Yang-Mills theory fl34p to 
compare the real-time evolution in the gravity. Because of the appearance 
of the cascade decay in the thermodynamics of the one dimensional model 
([1]), we expect that a similar decay pattern might appear in the real-time 
evolution, too. One aim of this calculation is to see the fractal structure 
during the evolution, which was found in the black string decay [IH]. This 
fractal structure is akin to a low viscosity fluid, and may indicate the ap- 
pearance of the naked singularity. Hence it plays a key role in the breaking 
of the cosmic censorship hypothesis. If we can observe the corresponding 
fractal structure in the gauge theory, it may be valuable in the context of 
the fluid/gravity correspondence [ID] and in the study of the resolution of 
the naked singularity through the quantum (finite- iV) effects [22l |33| IMf 135] . 
(In the stochastic time evolution of the model ([1]), any clear fractal structure 
has not been observed. We need further numerical calculations at larger 
to conclude it.) 

The last challenge is the numerical analysis of the decay process of the 
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smeared black DO brane solution (155]) in the IIA supergravity. The authors in 
Refs. [TSl dSl [2D] studied the decay of a 5 dimensional neutral black string in 
an asymptotically fiat spacetime with S^, which is distinct from the smeared 
DO brane geometry. To compare the real-time evolution in the gauge theory, 
it is important to study the DO brane geometry. 
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A General Multi-cut Solutions in the One- 
dimensional Gauge Theory 

In section 13. 2[ we derive Zm solution in the model ([T]). In this appendix, 
we consider the derivation of general multi-cut solutions at high temperature 



To derive an m-cut solution, we consider m mobs of the eigenvalues, where 
the l-th mob consists of A'^; eigenvalues. Here Ya^i = ^ assume 
Ni ~ 0{N) and m <^ N. Then we decompose the eigenvalues in the Z-th 
mob as 



T > T,2. 

The effective action f lTTj) can be rewritten as 




= - ^ . + - ^ log [2 cosh(/3A) - 2 cos( 
8A 2 . . 



«i - aj)] 




(39) 




ai = 9i + a) 



,(0 
j ' 



for J2Nk<i<J2^k, j=z-^A^, 



fc=l k=l k=l 



where 6i and represent the centre of mass and the relative position of the 
l-th mob. Oj"^ satisfies the constraint Ylf=i '^f^ ~ 0' which is negligible at 
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large N. We assume \af^\ ^ \6k — 6n\ and \af^\ ^ /3A ^ 1, and decompose 
the action (l39l) as 



m 



S'cfr(ai, A) = + + ^ , with (40) 

1=1 

Here •* and 6'/ appear only in 5*^(0 and Sq respectively, and we can solve 
them independently. 

First, we consider SaW- Then we see that af^ behave as a free Hermitian 
matrix model with a mass M = y/2D/l3A [27]. Thus af^ compose a cut with 
a width A/3A/y/2D and the assumption <^ /3A is justified at large D. 

Next we turn to S'^. m equations of motion are derived from Sg, 

f:^<=ot(^)=0. (/ = l.....m). (41) 

Although these equations are complicated, we would obtain a solution {Oi}, 
which satisfies 9k 9i, for a given set of Ni at least numericall\l^. 

Lastly we consider Sa- (Additional potentials for A arise by integrating 
af ■* in Sail) but they are negligible at large D.) has a unique saddle point 

-(|E(5)T" 

This result implies that the assumptions /3A ^ 1 and \af^\ ^ \dk — 9n\ are 
valid at sufficiently high temperature. In this way, we obtain the multi-cut 
solutions in the model ([T]). 



""^^ As far as we calculated, we always found at least one numerical solution of eq. (j41]l for 
a given set of Ni . 
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Lastly we evaluate the observables in this solution. The value of the 
effective action is given by 



pom f3Dm 





iV^Y log (/3A) 
Nj (3 



+ ■■■ 



(43) 



at the high temperature. The eigenvalue density becomes 



a 



(44) 



1=1 



1 

' i=i 



a 




One interesting feature of this density is that the widths of each cut are com- 
monly given by 4/M. Indeed, in the numerical calculation of the stochastic 
evolution in Fig. [TJ we observe that the widths of the cuts of steady states 
seem common. 

Note that A and 5*^ are determined by A''; and f3, and do not depend on 
9i. It may explain the reason that, in the stochastic evolution in Fig. [9l A 
remains almost constant until around s = 2200, whereas {u„} start varying 
from around s = 1000. In this evolution, the two mobs gradually approach 
each other, and finally they merge into one mob. Since {«„} are sensitive to 
the angle of these two mobs, they are varying during the evolution. On the 
other hand, since the numbers of the eigenvalues in each mob do not change 
until they merge, A may remain constant. 

B Details of the Monte Carlo Calculation 

We explain the details of the stochastic time evolution through the Monte 
Carlo calculation of the model ([T]), which we employ in Sec. HI 

Action: We adopt the same action introduced in Appendix B of Ref. [35] . 

Random number r: r in "1." in Sec. 14.11 is taken r = rA(0, 1), where 
r is a small number and A(0, 1) is a random number with standard normal 
distribution. Typically, we take r = 0.01. 
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Adjoint scalar: In order to achieve the equilibrium in the adjoint scalar 
at each step as argued "4." in Sec. 14. 1^ we repeat the heatbath updating 
of the scalar fields sufficiently many times. Here we observe that, if the 
repetition number is not enough, the stochastic evolution does depend on the 
lattice spacing. Smaller lattice spacing requires more repetitions to achieve 
the equilibrium. Typically we take the lattice spacing 0.05 and repeat 5 times 
at each step. About this lattice spacing, see Ref. [H]. We have confirmed 
that this repetition number is sufficient at this lattice spacing. 
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